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ABSTRACT 

We present the results of our analysis of six orbits of i^ST/ST/^ time-tagged 
spectroscopy of the intermediate polar, EX Hydrae. The high time and wave- 
length resolution of the HST/STIS spectra provided an excellent opportunity 
to study the UV properties of EX Hya. Measurements of the continuum and 
emission line fluxes corroborate earlier studies that show that the emission line 
fluxes are modulated more strongly than the continuum flux and originate from 
the accretion curtains, while the continuum flux originates from the white dwarf 
photosphere. The measured K amplitude of the narrow emission line radial ve- 
locity curve is used to calculate a white dwarf mass of 0.91 ± 0.05 Mq. Synthetic 
white dwarf photosphere and accretion disk spectral models are used to further 
refine the white dwarf and accretion disk properties. Using the spectral models, 
it is determined that EX Hya has a white dwarf of mass 0.9 Mq, T = 23000 K, 
an accretion disk truncated at 2.5 -Rwd, and is at a distance of 60 pc. 
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1. INTRODUCTION 

Intermediate polars (IPs) are a class of magnetic cataclysmic variables in which a white 
dwarf primary star sustains a magnetic field that is sufficient enough to divert the gravita- 
tionally controlled mass transfer. IPs may have a truncated accretion disk (the inner edge 
being wiped out by the magnetic field) or may be diskless accretors (accreting directly from 
the baUistic accretion stream), but the defining characteristic that separates IPs from the 
higher magnetic field polars is an asynchronously rotating white dwarf. As the magnetic 
field lines of the white dwarf sweep around the inner edge of the accretion disk, curtains 
of magnetically controlled material form and follow along the field lines until the material 
impacts the white dwarf surface at or near the magnetic poles, with both upper and lower 
poles in IPs accreting equally. 

EX Hydrae is an IP with a spin period of 67 minutes and an orbital period of 98 minutes, 
which causes an unusual scenario in which the corotation radius is close to the Li point (King 
& Wynn 1999). Magnetically controlled accretion occurs from the inner edge of the disk, 
and possibly from the outer edge of the disk, as evidenced by dips seen during the bulge 
dip in EUV light curves of EX Hya (Belle et al. 2002). Radiation from EX Hya is seen to 
be modulated on both the binary orbital and white dwarf spin periods. Binary modulation 
is seen in the form of a large bulge dip at phases (pgg — 0.5 — 1.1 (^gg refers to the binary 
phase) in the EUV light curve, which is caused by absorption along the line of sight to the 
emitting poles due to the extended hot spot on the edge of the accretion disk. Continuum 
emission is modulated sinusoidally over spin phase for all wavelength regimes (e.g., Hellier 
et al. 1987; Rosen et al. 1991; Belle et al. 2002). Maximum spin phase occurs when the 
upper accretion pole points away from the observer, allowing a direct view of upper and 
lower accretion regions; minimum spin fiux occurs when the upper pole points towards the 
observer and the lower pole is self-eclipsed by the white dwarf. 

Previous UV observations of EX Hya consist of HUT spectra (Greeley et al. 1997), 
ORFEUS II and /[/E spectra (Mauche 1999), and HST/FOS spectra (Eisenbart et al. 2002). 
These observations show that the broad emission lines and the continuum fiux are formed 
in different regions; the emission lines from an optically thick accretion curtain and the 
continuum from the inner most part of the accretion curtain and the white dwarf itself. 
Blue-shifted absorption was also seen in the C IV A1549A emission line. The UV data of EX 
Hya have also been modeled using a variety of continuum emission sources. The continuum 
has been fit by a T = 25000 K blackbody (Greeley et al. 1997), but this same work also 
showed that a bremsstrahlung model, with T = 74000 K, or a white dwarf photosphere 
of T = 26200 K also fit the data well. Eisenbart et al. (2002) found that a white dwarf 
photosphere of T^s = 25000 ± 3000 K (determined from line ratios) fit their FOS spectra 
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well. 



2. OBSERVATIONS 

Our HST/STIS observations, obtained through DDT time, comprised six orbits and 
used the E140M grating in Echelle mode. Each orbit produced a time-tagged spectrum 
consisting of 40 Echelle orders. The E140M grating spans a wavelength range of 1150 — 1735A 
and has a resolving power of = 45, 800, which corresponds to a wavelength resolution of 
0.013 — 0.019A/pixcl throughout the spectrum. The six orbit observation provided complete 
coverage of EX Hya's spin orbit, but only ~ 2/3 of its binary orbit. Table 1 gives the specific 
phase coverage of each data set. Each spectrum was processed with IRAF/STSDAS software 
and the time-tagged data was utilized to bin the spectra into phase resolved components. 

Spectra were extracted for a variety of orbital and spin phase bins. Spin phase bins 
were determined according to the spin modulation seen in the EUV light curve of EX 
Hya (Belle et al. 2002) and we have defined spin maximum phases as 06? = 0.865 — 0.365 
and spin minimum phases as <;/)(i7 — 0.365 — 0.865, where (f)^^ refers to the 67 minute spin 
ephemeris, T = 2437699.8914(5) + 0.046546504(9)E - 7.9(4) x lO^^^E^ (Hellier & Sproats 
1992). Spectra extracted throughout the spin period were binned into A^g? = 0.1 = 
402.2 s phase bins; bins were chosen as such to preserve a moderate signal-to-noise ra- 
tio. Three types of binary phased spectra were created: binary phased, spin maximum 
binary phased, and spin minimum binary phased. For the binary ephemeris, we used 
T = 2437699.94179+0.068233846(4)E (Hellier & Sproats 1992). Spin maximum data existed 
for binary phases (^gg = 0.0 — 0.7, but spin minimum data only existed for binary phases 
098 = 0.0 — 0.6. Each binary phase bin has a resolution of A^g^ = 0.1 = 589.5 s, where 
098 refers to the 98 minute orbital period. We used IRAF/SPLOT to analyze the spectra 
and determine the la errors we report here. Figure 1 presents our complete spectral data 
set divided into two phase bins: spin maximum and spin minimum. The spin maximum 
spectrum has on overall higher continuum level and stronger emission lines than the spin 
minimum spectrum. Each spectrum exhibits broad emission lines of C III A1175A, Si III 
A1206A, A1298A, Lya A1216A, N V A1240A, C II A1335A, O V A137lA, Si IV A1396A, C IV 
A1549A, and He II A1640A, along with narrow emission lines N V AA1238, 1242A and O V 
A1371A, which appear in all phases except for the early binary phases (09g = 0.0 — 0.1). 



-4- 



3. DISCUSSION AND ANALYSIS 
3.1. Continuum 

Continuum flux measurements were made at three emission line-free wavelength regions 
along each spectrum: AA1352 - 1362A, AA1440 - 1520A, and AA1565 - 1620A. The results 
of the AA1352 — 1362A measurements are shown in Figure 2 (spin phased). The other two 
wavelength regions exhibit similar behavior. Error bars on each measurement are RMS values 
as reported by IRAF/SPLOT. The spin phased continuum flux (Figure 2) is fit well by the 
sinusoidal function, /(10~^^ ergss""*^ cm^^ A^-*^) = A + i?sin27r(067 — 0o)) where A = 1.68, 
B = 0.38, and (f)o = 0.76 ±0.05, which places the peak of the light curve at = 0.01 ±0.05. 
The error here is dominated by the fact that each phase bin measurement was confined 
within the A^g? = 0.1 phase bin resolution of the data. The continuum flux is modulated by 
B/A — 23%. The peak in UV continuum flux is earlier than that of the EUV emission, which 
occurs at = 0.115 ± 0.001 (Belle et al. 2002). As both the HST and EUVE observations 
of EX Hya occurred simultaneously, an error in the ephemeris is ruled out as a cause for 
phase difference in peak continuum fluxes; rather, the phase disparity is likely a result of 
distinct EUV and UV continuum emission sites. The sinusoidal fits to all of the spin phased 
UV continuum wavelength regions are given in Table 2. 

The spectra were also measured over binary phase, and, in agreement with previous 
studies (e.g., Hellier et al. 1987), wc found that the continuum flux did not exhibit coherent 
modulation over the binary orbital period. 

3.2. Emission Lines 

3.2.1. Broad Emission Lines 

The broad emission lines were fit with single Gaussian functions and a linear continuum, 
using IRAF/SPLOT, to determine parameters such as central wavelength, flux, FWHM, and 
equivalent width. The Lya and Si III A1206A lines were not measured because they are 
blended during the spin maximum phases. 

Over spin phase, the emission line fluxes exhibit a sinusoidal shape (C IV is shown 
as an example in Figure 3), and can be flt well with a sinusoidal function of the form 
/(10~^^ergss~^ cm~^) — A + 5sin27r(067 — 0o)- The solutions for each line are given in 
Table 3. The emission peaks near 067 ~ 0.05 — 0.1 for each line. Considering the error of 
0.05067 on each phase bin, the emission line fluxes peak at roughly the same phase as the 
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continuum flux (067 ~ 0). C III A1175A has the greatest modulation, at 73%, but aU of 
the hues are strongly modulated by around 50%. Compare this with the continuum with 
a flux modulation of only ~ 20 — 25%. While both the continuum and emission lines are 
modulated sinusoidally over spin phase, the difference in the modulation amplitude suggests 
that the continuum and emission lines have different formation sites. Over binary phase the 
emission line ffux is not modulated coherently. Comparison of the emission line fluxes over 
binary and spin phase demonstrates that the flux modulations are clearly a function of spin 
phase, rather than binary phase. 

Equivalent widths (EWs) of strong UV emission lines from optically thin plasmas are 
typically on order of 100s of Angstroms. The EWs of the emission lines appearing in the EX 
Hya spectra are smaller than expected for optically thin plasma emission. C IV, which has 
the largest EW of all the broad emission lines, is shown as an example in Figure 4 (solutions 
for the sinusoidal flts to all of the hnes are given in Table 4) . There are two possible scenarios 
that may explain the observed EWs. The first is that optically thick material is responsible 
for the emission lines. The second is that the optically thin material is only weakly emitting 
and another form of radiation, such as cyclotron emission, may play an important role in the 
radiative cooling of the accretion curtains. If the latter is true, evidence should be seen in 
the form of cyclotron humps; for a magnetic field strength of 5 ~ 5 MG, the first cyclotron 
harmonic should appear around 10/x. The almost triangular shapes of the lines implies that 
they are not from a purely optically thin material, thus it is likely that both optically thin 
and optically thick emission contributes to each emission line. 

Figures 5 and 6 present two examples of the broad emission line evolution over binary 
and spin orbit. The N V A1240A fine shown in Figure 5 has been separated into spin 
maximum (thick line) and spin minimum (thin line) and phased over the binary orbit. The 
line retains roughly the same form throughout the binary orbit for both spin phases. Aside 
from binary phase bin A^gg = 0.0 — 0.1, the flux of the spin maximum line is constantly 
greater than the flux of the spin minimum line. This suggests that the lines are formed in 
a region that would be partially absorbed or occulted during spin minimum phases, i.e. the 
inner part of the accretion curtain. Also of interest is the lack of narrow emission lines during 
binary phases A^gg = 0.0 — 0.1, which correspond to the egress of the bulge dip; absorption 
features appear in their place. The narrow lines will be discussed in further detail in §3.2.2. 
The fact that the hnes have roughly the same amount of absorption during A^gg = 0.0 — 0.1 
implies that the lines are formed interior to the bulge. 

If the change in flux of the broad emission lines between spin phases maximum and 
minimum is due to absorption through the accretion curtain, we can estimate the density 
of the curtains from the amount of absorption. Using the simple relation of / = Io^~^, 
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we may first deduce r and then assume a Kramer's opacity law to determine the density. 
The physical size of the accretion curtains must first be estimated. At most, the radial 
extent of the footprint of the curtain on the disk could be the entire truncated accretion 
disk and at the least, it would only be a small fraction of the disk size. The disk size was 
determined by measuring optical emission lines (using the optical data obtained by us during 
this campaign) to determine the inner and outer accretion disk radius, assuming a Keplerian 
accretion disk. We considered two white dwarf masses (see §3.3), My/o ~ 0.90 M0 and 
MwD — 1.33 Mq, and calculated the disk radii corresponding to each white dwarf mass 
(-Rinner = 1.6 — 2. 3x1 0^ cm and -Router = 2.5 — 5. 9x1 0^ cm). Several different physical sizes, 
of the accretion curtain were used for our density calculation. We take 5" = o-Rdisk) where 
a = 0.01, 0.10, 0.50 and 1.0. In this way we were able to sample a range of physical depths 
for the accretion curtain. The emission line fluxes were measured during spin maximum 
and spin minimum, excluding binary phases of the bulge dip. This ensures that material 
associated with the bulge, which is also responsible for absorption at certain spin phases, 
will not be mistakenly identified with accretion curtain material. Five emission lines were 
used for this calculation; C III, N V, Si IV, C IV, and He II. Si III A1296A and C II were 
not used due to their ambiguous line profiles during spin minimum. For each white dwarf 
mass and corresponding S values, we find that r = 0.7 — 1.1 and the density ranges between 
XQ-io _ io~^ g cm^"' for all values of 5*. If we assume that the absorbing material is pure 
hydrogen, we obtain a number density of A^h = 7 x 10^^ — 8 x 10^"^ cm~^ within the accretion 
curtains. 

The C IV A1549A emission line has a blue-shifted absorption component that appears 
during spin phases A^gy ~ 0.2 — 0.8 as evidenced in Figure 6. A similar feature was also 
observed in the C IV line in HUT data during spin phases A^g^ = 0.36 — 0.48 (their data 
only extended through = 0.5, Greeley ct al. 1997). Although their C IV line is phased 
into finer orbital bins, our C IV line clearly shows an absorption component that begins prior 
to 067 = 0.3 and evolves into an absorption that eats away the blue wing of the emission 
line. At phases 067 — 0.2 — 0.3, the blue absorption has a FWHM of ~ 275kms~^ and is 
blue-shifted by ~ 615kms~^ from the center of the C IV doublet. This is slightly less than 
the value of ~ 1000 km s~^ reported by Greeley et al. (1997) for their blue G IV absorption 
feature. 

G III, N V, and G II also show a blue absorption component during spin phases 067 ~ 
0.2 - 0.7. The N V absorption component has a FWHM of ~ 400kms~^ and is blue-shifted 
by ~ 700 km s"^ from the center of the line. The G III absorption component is really too 
shallow to accurately measure, though it is visible during spin phases 067 ~ 0.2 — 0.7. The 
G II line seems to have a blue-shifted absorption component throughout the entire spin phase. 
Its FWHM ranges from 655 km s"^ at 067 = 0.4-0.5 to 200 km s"^ at 067 = 0.0-0.1 while 
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the blue-shift from the center of the emission hne remains fairly constant at ~ 1400 km s~^. 
However, the C II absorption feature is associated with a white dwarf absorption line (see 
§3.4), which may influence the motion of the line. 

Radial velocity curves for the broad emission lines are shown in Figure 7 phased over the 
spin period. The spin phased radial velocity curves present surprising results. In previous EX 
Hya observations (e.g., HcUicr ct al. 1987), the maximum blue-shift of the optical emission 
lines observed around (j)i„ = 0.0 has been used to place the origin of these lines in the 
accretion curtain; while the upper pole points away from us, material flowing into the pole 
will be at its maximum blue-shift. We sec that three lines exhibit this behavior: C III, N V, 
and C II. Si III, Si IV, and He II, on the other hand, have maximum blue-shifts closer to 
067 ~ 0.5. A double peaked structure is visible in all but the Si III radial velocity curve. 



3.2.2. Narrow Emission Lines 

Three narrow emission lines appear in the HST spectra: N V AA1238.8, 1242.8A and 
O V AI37I.3A, and represent the two highest ionization energy lines present in the spectrum, 
with XNeV = 98 eV and xov = 114 eV. These lines have FWHM~ O.2A and are visible over 
both spin and binary orbit, and only disappear during binary phases = 0.0 — 0.1, which 
correspond to the egress of the bulge dip. A radial velocity curve over spin phase is shown 
for the narrow emission lines in Figure 8. Each radial velocity curve has a double-peaked 
structure, although the O V and N V A1238A lines are less modulated than the N V A1242A 
line. For each line the velocities tend to be blue for two reasons: one, only two-thirds of 
the binary orbit, (pgg = 0.0 — 0.7, are (unevenly) sampled in each spin phase bin and two, a 
7— velocity has not been subtracted. Previous measurements of the 7— velocity (derived from 
optical spectroscopy) range from — 9±35kms~^ (Gilliland 1982) to — 180kms~^ (Hellier et 
al. 1987), and we will show below that we obtain a value shghtly different from the GiUiland 
(1982) and HeUier et al. (1987) results. 

The radial velocities of the narrow emission lines over the binary phase arc shown in 
Figure 9 and can be fit with a sinusoidal function of the form v = j + sin27r(098 — (po), 
with 7 = 9.5 ± 3kms-^ K = 59.6 ± 2.6 kms"^ and 00 = 0.98 ± 0.05. The K amplitude 
derived here is comparable with a previously determined value of Ki — 69±9kms~^ (Hellier 
et al. 1987). The small velocities, the clear sinusoidal modulation of the velocities, and the 
K amplitude suggest that these lines are formed near or on the white dwarf surface. A high 
temperature emission region for these lines increases the likelihood that the lines are not 
contaminated by emission from the accretion disk. This narrow line radial velocity curve is 
therefore likely to be a true representation of the motion of the white dwarf, as compared to 
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optical (disk) emission lines that have been used previously to determine the Ki amplitude 
and the 7— velocity (e.g., Hellier et al. 1987). We would like to note that the proximity of 
the narrow line emission region to the white dwarf surface should cause these lines to be 
gravitationally redshifted. Unfortunately, without an accurate 7-velocity measurement, it is 
impossible to separate the gravitationally redshift from the 7-velocity of the narrow emission 
line radial velocity curve. 

The narrow emission lines appear to be double-peaked during certain spin phases. The 
N V doublet shows a strong double-peaked profile during spin phases 06? = 0.1 — 0.3, 
067 = 0.5 — 0.6, and ({)(„ = 0.9 — 1.0. O V is double-peaked during spin phases (/)(„ = 0.0 — 0.3, 
067 = 0.4 — 0.6 and 06? = 0.7 — 0.8. These phases roughly correspond to spin maximum 
phases. Close inspection of Figure 5 also shows the double-peaked nature of the N V line 
during spin maximum and that during spin minimum, it is only the blue component of the 
double-peaked hne that is still visible. The radial velocities of these hues over binary phase 
tell us that they originate close to the white dwarf surface. The double-peaked nature of 
the lines suggest a two velocity component as their origin. We propose that these lines are 
formed in a wind or material outflowing from the poles. During spin maximum, we have a 
direct view of both poles (possibly not a direct view of the northern pole as it points away 
from us, but a wind would clearly be seen) giving us both red and blue-shifted components. 
At spin minimum, when the upper pole points towards the observer, the lower pole and 
its (red-shifted) emission is blocked from view; only the blue-shifted narrow line component 
would be seen at these phases. As discussed in §3.2.1, evidence of a wind is also seen in the 
form of blue-shifted absorption in several of the broad emission lines. 



3.3. Mass Determinations 

In order to calculate a white dwarf mass, a value for the secondary mass must be 
estimated. Previous studies have used an M2 value determined from a zero age main sequence 
(ZAMS) mass-radius relation apphed to the Roche lobe radius for a given binary orbital 
period. HeUier et al. (1987) found M2 = 0.128 M© using the ZAMS mas^radius relation 
given in Patterson (1984) and an earlier study used the relation presented in Warner (1976) 
to find M2 = 0.18 M0. We use a value based on an up-to-date mass-radius relation (Howell 
et al. 2001) expressed in the form of 

M, = omr'-^'p^i (1) 

where Porb is the binary orbital period in hours and / is a factor that describes the bloating 
of the secondary star in systems above the period gap. For the EX Hya orbital period of 1.6 
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hr and f — 1, we calculate M2 — 0.152 M©, a value slightly higher than the 0.128 that 
has been used previously to determine the white dwarf mass. 

Using the Ki amplitude determined from our UV narrow emission line radial velocity 
curve, we expect to obtain an unambiguous determination of the mass of the white dwarf if 
our value of M2 is correct. We employ a modified form of Kepler's third law, 



and solve for q in order to determine Mw^. For i = 78 ± 1° (Hellier et al. 1987), M2 = 
0.152 Mq, and Ki — 59.6 ± 2.6kms~^, we determine a white dwarf mass, Mwd — 1-33 ± 
0.11 Mq. This value is considerably higher than 0.78 ± 0.17 M©, which was determined by 
HeUier et al. (1987), and the even lower values of ~ 0.5 Mq that have been determined from 
X-ray studies (e.g., Fujimoto & Ishida 1997). Although our Ki amplitude agrees within the 
errors with the value measured by Hellier et al. (1987), our white dwarf mass does not agree 
with theirs due to the larger secondary mass used in the above calculation. 

Smith, Cameron, & Tucknott (1993) apphed the method of skew-mapping to detect the 
secondary stars in cataclysmic variables. Applying this to EX Hya, they found a velocity 
amplitude for the secondary motion of K2 = 356 ±4kms~^. This amplitude, along with our 
measured value of Ki and calculated value of M2 may be used to apply another determination 
of Ml using K1/K2 = M2/M1. This method gives a white dwarf mass of Myy^ = 0.91 ± 
0.05 Mq, which is ~ 30% lower than the value determined above, although still larger than 
the earher estimates (e.g., HeUier et al. 1987). 



Discerning exactly which components of the binary system are contributing UV contin- 
uum flux (and at which spin phases) is an important issue for thorough modeling of the HST 
spectra of EX Hya. Spin maximum phase would seem to be the most likely time in which 
we would have a direct view of the white dwarf, but are there other components also visible 
during this phase? The accretion poles (shock regions) are strong emitters in the X-ray 
and EUV (Belle et al. 2002). The irradiated material in the accretion curtain close to the 
pole, however, may emit reprocessed X-ray and EUV radiation as UV radiation, thus adding 
another UV emission component to the system. The accretion disk, although truncated, is 
likely to add to the UV continuum emission as well. 

To model our HST spectra with a white dwarf photosphere, we must determine which 
spectrum is the best representation of the white dwarf by ascertaining which of the different 
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3.4. Spectral Modeling 
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components discussed above contribute to the spin maximum and spin minimum spectra. 
Visual inspection of the spin maximum and minimum spectra reveals their differences; the 
red- ward slopes of each spectrum are roughly the same while the blue end of the spin min- 
imum spectrum has a lower flux than that of the spin maximum spectrum. There are two 
possible cases to consider; 1— spin minimum is the true white dwarf spectrum and the ac- 
cretion curtain is optically thin to the white dwarf continuum. Spin maximum therefore 
contains an additional continuum component that could be the accretion region located at 
either pole. 2— The spin minimum spectrum is absorbed by the accretion curtains and the 
spin maximum spectrum is the true white dwarf spectrum plus an additional continuum 
component. These two options outhne the occultation (1) versus absorption (2) models. For 
option 1, spin phased modulation is due to occultation of an emitting pole, while option 
2 dictates that spin phased modulation is due to UV flux being absorbed by the accretion 
curtains. 

For either model option, the spin maximum spectrum has added UV continuum compo- 
nents. This is evidenced by the fact that the blue slopes of the spin maximum and minimum 
spectra are not equal. If the spin modulation were simply due to absorption, the spin mini- 
mum spectrum would retain the same slope as that of the spin maximum spectrum. (This 
assumes a gray absorber, which is valid over the small, ~ 600A, wavelength range covered by 
the //ST spectra.) We show the difference spectrum of EX Hya, which is the spin maximum 
minus the spin minimum spectrum, in Figure 10. This spectrum represents the additional 
UV continuum and emission line components that are seen during spin maximum phases. 
Considering the above arguments, we rule out the spin maximum spectrum as representa- 
tive of only the UV white dwarf continuum flux. For option 1 to be viable, the accretion 
curtains near the poles would have to be dense enough to absorb the flux emitted from the 
accretion region. While the outer curtains arc optically thin, it is likely that the curtain 
close to the pole is optically thick and will absorb emitted radiation from the accretion pole. 
We therefore suggest the case that the spin minimum spectrum is the better representation 
of the pure white dwarf spectrum and that a blue blackbody and emission line component 
(or multiple components) is seen in addition to the white dwarf during spin maximum. Flux 
contributed by the truncated accretion disk must also be considered at all spin phases. 

We first created white dwarf photosphere models using TLUSTY and SYNSPEC (Hubeny 
1988) for the two white dwarf masses computed in §3.3, 0.90 (log g = 8.5) and 1.33 
(log g — 9.4). We modeled white dwarf photospheres of solar composition and considered 
temperatures ranging from T^s — 19000 K — 35000 K (in steps of 1000 K). Although it was 
decided that the spin minimum spectrum would be the best representation of the white dwarf 
spectrum, the model photospheres were fit to the spin maximum and minimum spectra in 
order to insure that the spin minimum spectrum was indeed the best choice. For each fit. 
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the broad emission lines in the spectrum were masked, except for the wavelength range 
including the Si IIIA1298A and C IIA1335A emission lines because of the white dwarf ab- 
sorption components that appear between these emission features. Our best fitting 0.90 Mq, 
Tgff = 23000 K white dwarf photosphere model is shown in Figure 11 fit to both the spin 
maximum and spin minimum spectra. The spin minimum fit has a reduced of 1.7, and 
scaled to the i/^T continuum fiuxes, gives a distance to EX Hya of 54 pc. The spin maximum 
fit has a reduced value of 3.3 and scaled to the EST continuum fluxes gives a distance of 
46 pc. 

Note that many of the absorption features appearing between the Si III and C II in the 
spin minimum spectrum are fit well by the white dwarf photosphere model. The slope of the 
model spectrum also produces a better match to that of the spin minimum spectrum than 
the maximum spectrum. These facts, coupled with the lower reduced value of the spin 
minimum fit, corroborate with our earlier statement that the spin minimum phases are the 
ideal phases during which to model the white dwarf photosphere. We therefore use only the 
spin minimum spectrum in our following model fits. 

White dwarf photospheres with temperatures of Teflf = 23000 — 35000 K (in steps of 2000 
K) were created for M^^d = 1.33 Mq. The best fit model to the spin minimum spectrum has 
a temperature of Tcfj = 25000 K, but while this model provides a good fit, it also places the 
distance of EX Hya at 30 pc! This distance varies greatly from the distances of 65 — 100 pc 
that have been determined in previous studies (e.g., Hellier et al. 1987; Eisenbart et al. 2002). 
We explored the relationship of distance versus white dwarf photosphere temperature for 
both white dwarf masses over a range of temperatures and found that, for Mvkd = 1.33 M0, 
a white dwarf temperature of Teg = 35000 K would be required to achieve a distance of 55 pc 
(when scaled to the spin minimum spectrum). An My^D — 1.33 Mq, Teg = 35000 K white 
dwarf photosphere model was considered for the spin minimum spectrum, but the slope did 
not match the iJST continuum slope. However, before we completely rule out either white 
dwarf mass, we must also consider the accretion disk contribution to the UV continuum fiux. 

Accretion disk models were created using TLUSDISK, SYNSPEC, and DISKSYN (Hubeny 
1988; Wade h Hubeny 1998) for white dwarf masses of M^^d = 0.90 Mq and Mv^d = 
1.33 Mq. We used an accretion rate of lO^^gs"^ (Fujimoto & Ishida 1997) and an inclination 
of i = 78° (Hellier et al. 1987). In order to determine the radius at which the disk should be 
truncated, the inner disk radius was determined from our optical spectra. For a Keplerian 
disk, the highest velocities in an emission line originate from the inner disk. The optical line 
wings, therefore, should map to the inner edge of the disk. The Ha A6563A line profile in 
each of our optical spectra was used to determine the velocities at the blue and red wings. 
For a MwD — 0.90 white dwarf, these velocities map to an inner accretion disk radius 
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of -Rinner = 2.5 i?wD (1-6 X 10^ CHi) and for a Mwd — 1.33 white dwarf, the inner radius 
becomes i^inner = 9.0 /?wD (2.3 X 10^ cm). Previous modehng studies (e.g., Greeley et al. 
1997) have assumed -Rinncr ~ 5 i?wD (3.5 x 10^ cm) and optical spectroscopy by Hellier et al. 
(1987) gave i?inner = 6 X 10^ cm for Mwd ~ 0.8 M©. 

Two different accretion disk models must be created for the two different white dwarf 
masses we have been using. M\yd = 0.90 Mq and M = lO^^gs"^ correspond to Models 
n and p of Wade & Hubeny (1998) and M^d = 1-33 Mq with the same accretion rate 
matches most closely with Models y and z. Model n {p) has UV flux contributing annuli 
that extend to 6.2 -Rwd (9.40 Rwd), and Model y (z) has contributing annuli that extend out 
to 14.3 -Rwd (21.7 iiwo); the last annulus in each model has T ~ 10000 K. Although it would 
seem that the synthetic spectra for each disk model would be different because of the different 
truncation points, the truncated accretion disk spectra are in fact quite similar. This is due 
to the fact that the more massive white dwarf has a hotter accretion disk and while the disk 
is truncated at a greater radius than the less massive white dwarf, an approximately equal 
amount of UV flux is still emitted from the remaining annuli. 

The truncated disk models were added to the white dwarf photosphere models using 
the following method: the photosphere model and disk model were scaled separately to a 
distance, d, and then summed together. The summed spectrum was then fit to the observed 
spectrum. The distance, d, was modified until the summed spectrum produced a good fit 
to the observed spectrum. The summed model fits were judged by visual inspection, rather 
than by a formal fitting routine. For Myy^ = 0.90 Mq, we found that a T^s = 23000 K 
white dwarf plus the truncated disk produced the best fit at a distance d = 60 pc. The 
Mwd = 1.33 M0 model showed a best fit for T^s = 25000 K plus the truncated disk at a 
distance of d = 33 pc. The M-wd — 0.90 Mq fit is shown in Figure 12, as we definitively rule 
out the higher mass model below. 

The white dwarf photosphere plus truncated accretion disk model fits to the HST fiuxes 
give two different distances: the 0.90 Mq white dwarf model gives d — 60 pc, while the 
1.33 Mq white dwarf gives d — 33 pc. We can check these distances independently by using 
Bailey's method (Bailey 1982), which uses the band magnitude and measured surface 
brightness values for K-M single stars to determine a distance; 



where Sk is the surface brightness. Using K — 11.705(28) from the 2MASS survey (Hoard et 
al. 2002), Sk = 4.5±0.5 for an M3 main sequence star (Dhillon et al. 1997), R2 = 0.17i?Q, and 
the secondary band flux contribution as 75 ±25% (Dhillon et al. 1997), we determine the 
distance to EX Hya to be d = pc. This effectively rules out the 1.33 Mq, T^s = 25000 




(3) 



-13- 



K white dwarf, as the distance determined by this model is 33 pc. 

The greatest uncertainties in the distance calculation and the white dwarf mass deter- 
mination are the secondary star mass and spectral type. An M3 main sequence secondary 
star (determined from IR spectra, Dhillon et al. 1997) has a mass of ~ 0.3 M©. This is 
not consistent with the mass determined using the ZAMS mass-period relation (Equation 
1), which gave the secondary mass as 0.15 Af0. Returning to Equation 2, we find that a 
secondary mass of 0.12 Mq would be required to obtain M^^o = 0.9 Mq. A more accurate 
secondary mass will be needed to aid future distance and white dwarf mass determinations 
of EX Hya. 

4. CONCLUSION 

We have made an accurate measurement of the K amplitude of the white dwarf based 
on the narrow emission line radial velocity curve. The small velocities and high ionization 
temperatures required for production of the N V and O V emission lines indicate that they 
are formed close to the white dwarf, and therefore, the Ki = 59.6 ± 2.6kms~^ value that 
we have measured is an accurate representation of the motion of the white dwarf within the 
binary system. 

Synthetic spectral models for a white dwarf photosphere and truncated accretion disk 
were made using TLUSTY and SYNSPEC (Hubeny 1988). Based on spectral fits using a 
white dwarf photosphere model, we determined that the spin minimum phases are the best 
representation of the white dwarf spectrum, i.e., occultation of (one of) the emitting poles is 
responsible for the continuum modulation seen over the white dwarf spin phase. The distance 
we derive based on spectral modeling, d ~ 60 pc, agrees well with two recent measurements; 
Eisenbart ct al. (2002) determine c? = 65 ± 11 pc using Bailey's method and a K magnitude 
derived from their infrared spectra, while a recent parallax measurement has placed the 
distance to EX Hya at ~ 60 pc (K. Beucrmann, private communication). This implies 
that our spectral modeling parameters for a solar composition white dwarf photosphere and 
truncated accretion disk - Mwd = 0.9 M0, T = 23000 K, and -Rinner = 2.5i?wD - represent 
EX Hya well. 

We thank the director of the Hubble Space Telescope for allocating DDT time for the 
EX Hya multi- wavelength project. We also thank K. Beuermann for use of his EX Hya 
parallax distance prior to publication. The anonymous referee is thanked for very construc- 
tive comments. Support for proposal 8807 was provided by NASA through a grant from 
the Space Telescope Science Institute, which is operated by the Association of Universities 
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Table 1. HST spectra, phase coverage. 



HST Data set 


HJD start-sto]) 


S])in 


Phase 


Binary 


Pliase 


O68301010 


2451682.894 - 


2451682.919 


0.67 


-0.22 


0.32- 


■0.69 


O68301020 


2451682.950 - 


2451682.982 


0.88 


-0.55 


0.14- 


■0.60 


O68301030 


2451683.017- 


2451683.049 


0.32 


-0.99 


0.13- 


■0.59 


O68302010 


2451694.622 - 


2451694.647 


0.62 


-0.17 


0.19- 


■0.56 


O68302020 


2451694.680 - 


2451694.712 


0.87 


-0.55 


0.04- 


-0.50 


O68302030 


2151691.717- 


2151691.779 


0.32 


- 0.99 


0.03 - 


0.18 



Table 2. Sinusoidal fits to the continuum fluxes. 



Continuum Region 0o B/A 

AA1352 - 1362A 1.683 d= 0.001 0.381 ± 0.001 0.77 ±0.05 23% 

AA1440 - 1520A 1.554 ±0.001 0.315 ±0.001 0.76 ± 0.05 20% 

AA1565 - 1620A 1.375 ±0.001 0.337 ±0.001 0.78 ±0.05 24% 



'^A and B have units of 10 ergss ^ cm ^ A ^. 
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Tablc 3. Simisoidal fits to the broad emission line fiuxes. 



Ion 






00 


B/A 


cm 


1.897 ± 0.004 


1.381 ± 0.003 


0.81 ±0.05 


73% 


N V 


1.694 ±0.002 


0.864 ±0.001 


0.80 ±0.05 


51% 


Si III 


1.183 ±0.001 


0.618 ±0.001 


0.83 ±0.05 


52% 


CII 


1.288 ±0.001 


0.555 ±0.001 


0.81 ±0.05 


43% 


Si IV 


2.895 ±0.002 


1.732 ±0.001 


0.81 ±0.05 


60% 


CIV 


6.771 ±0.004 


3.490 ± 0.002 


0.80 ±0.05 


52% 


He II 


0.963 ± 0.001 


0.554 ±0.001 


0.82 ±0.05 


58% 



and B have units of 10 ergss ^ cm 



Table 4. Sinusoidal fits to the broad emission fine equivalent widths. 



Ion 


^a 


B^ 


00 


B/A 


C III 


10.84 ±0.02 


-5.35 ±0.01 


0.81 ±0.05 


49% 


N V 


10.44 ±0.01 


-1.738 ±0.001 


0.82 ±0.05 


17% 


Si III 


7.20 ±0.01 


-2.141 ±0.002 


0.90 ±0.05 


30% 


Si IV 


17.06 ±0.01 


-6.638 ± 0.003 


0.82 ± 0.05 


39% 


CIV 


47.03 ± 0.03 


-12.37 ±0.01 


0.83 ± 0.05 


26% 


He II 


7.39 ±0.01 


-2.078 ±0.002 


0.85 ± 0.05 


28% 



^A and B have units of A. 




Fig. 1. — Our HST/STIS data set separated into spin maximum and spin minimum phases. 
The spectra contain broad emission hues of C III A1175A, Si III A1206A, A1298A, Lya 
A1216A, N V A1240A, C II A1335A, O V A137lA, Si IV A1396A, C IV A1549A, and He II 
A1640A, along with narrow emission hues N V AA1238, 1242A and O V A137lA. The thin 
hne in both plots represents the Icr errors. 
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Spin Phase 



Fig. 2. — Continuum flux measured at AA1352 — 1362A plotted over spin phase and fit with 
a sinusoidal function of the form /(10~^^ ergss"^ cm~^ A^^) = A + B sin 27r(0B7 — (po), where 
A = 1.683 ± 0.001, B = 0.381 ± 0.001, and 0o = 0.76 ± 0.05. 



- 20 - 




Spin Phase 



Fig. 3. — C IV broad emission line flux plotted over spin phase and fit with a sinusoidal 
function (solution given in Table 3). 
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0.5 1 1.5 2 

Spin Phase 



Fig. 4. — Equivalent width of the C IV broad emission hne phased over the spin period. 
Solutions for the sinusoidal fits to all of the broad emission lines are given in Table 4. 
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1230 1235 1240 1245 1250 1230 1235 1240 1245 1250 1230 1235 1240 1245 1250 1255 



Wavelength (A) 



Fig. 5.— N V A1240A emission line profile evolution over the binary orbit and separated 
into spin maximum (thick line) and spin minimum (thin line). For all phases except 0gg = 
0.0 — 0.1, the line profile extracted at spin maximum is greater in flux than the spin minimum 
line. Also note the narrow emission lines that appear during phases = 0.1 — 0.7 and 
that exhibit a double-peaked nature in the spin maximum spectrum near binary phases 
098 = 0.1 — 0.2, while the blue component only is visible in the spin minimum line profile. 
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Fig. 6.— C IV A1549A line over spin phase. Note the blue-shifted absorption that appears 
during phases = 0.2 — 0.3 and the asymmetric profile as the line evolves over spin phase. 
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Fig. 7. — Broad emission line radial velocities phased on the spin period. Si III is the only 
line that appears to have a sinusoidal shape; the other radial velocity curves are somewhat 
double-peaked, with one peak occurring near = 0.2 and the second occurring at — 0.5. 
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Fig. 8. — Narrow emission line radial velocities phased on the spin orbit. The radial velocity 
curve of each line is double-peaked, with maxima occurring at = 0.2 and = 0.8; the 
N V A1242A component shows the strongest modulation. 
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Binary Phase 

Fig. 9. — Binary phased radial velocity curve for the narrow emission lines. This radial 
velocity curve is fit nicely with a sinusoidal function of the form w = 7 + ii'sin27r(0g8 — (Pq), 
with 7 = 9.5 ± 3kms-\ K = 59.6 ± 2.6 km s'^ and 0o = 0.98 ± 0.05. 
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Fig. 10. — Spin maximum - spin minimum difference spectrum. 
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Fig. 11.— An MwD = 0.90 Mq, T^s = 23000 K white dwarf photosphere model (black 
spectra) fit to the spin maximum (top) and spin minimum (bottom) spectra (gray spectra). 
The slope and absorption lines of the spin minimum spectrum are fit better by the white 
dwarf photosphere model than the spin maximum spectrum. 
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Fig. 12. — Best fit white dwarf photosphere plus truncated accretion disk model for EX Hya. 
This model uses Mwd = 0.90 Mq, T = 23000 K, and i?inner = 2.5 -Rwd, and is fit to the spin 
minimum spectrum. Scaled to the EX Hya spectra, this model gives a distance of 60 pc. 



